Absence of Anderson localization of light in a random ensemble of point scatterers 
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As discovered by Philip Anderson in 1958, strong disorder can block propagation of waves and 
lead to the localization of wave-like excitations in space. Anderson localization of light is particularly 
exciting in view of its possible applications for random lasing or quantum information processing. We 
show that in contradiction to the common belief, Anderson localization of light cannot be achieved 
in a random three-dimensional ensemble of point scattering centers that is a paradigm model to 
study multiple scattering of waves. This suggests that fully disordered systems are not suitable 
for observation of Anderson localization of light and casts doubt on the very possibility of this 
phenomenon for light in any three-dimensional medium without sufficient structural correlations. 



Anderson localization — the appearance and dominance 
of localized states in strongly disordered systems — is be- 
lieved to be a universal phenomenon for all quantum and 
classical waves p~H3]. In particular, three-dimensional 
(3D) disordered systems are expected to exhibit a tran- 
sition from the "metallic" phase with extended states to 
the "insulating" one with localized states, upon increas- 
ing the disorder. This transition was observed for elec- 
trons in disordered solids [4 j, ultrasound [5], and cold 
atoms [SHE] . Reports of Anderson localization of light in 
3D also exist [9tiTT] . Here we present a theoretical study 
of light scattering in a 3D ensemble of resonant point 
scatterers (atoms) at random positions. We show that 
Anderson localization takes place only in the scalar ap- 
proximation and disappears when polarization effects are 
taken into account. Resonant point scatterers allow for 
the strongest overall scattering thanks to the combina- 
tion of a large scattering crosssection and the possibility 
to pack an arbitrary large number of scatterers in a finite 
volume. Therefore, our results question the very possi- 
bility of Anderson localization of light in any 3D random 
medium, at least as far as there are no correlations in 
scatterer positions. 

Samples used in experiments on Anderson localization 
of light — semiconductor [9j [12] or dielectric [10] [TT] pow- 
ders, porous semiconductors [13] , or disordered photonic 
crystals [14] — are complex photonic media that do not 
allow for accurate theoretical description of light scatter- 
ing. Instead, theoretical physicists like to work with a 
simpler model, an ensemble of point scatterers, and im- 
plicitly assume that the overall transport properties of 
large disordered samples are independent of their micro- 
scopic details. Such an approach can indeed be quite 
successful and allows for a good understanding of un- 
derlying physics [15] [16]. In addition, the point-scatterer 
approximation is excellent for description of light scatter- 
ing on ensembles of cold atoms which therefore provide a 
fantastic and practically realizable playground for testing 
the theory [T7] . 

Let us apply the point-scatterer model to study Ander- 
son localization of light and try to go as far as possible 



without additional approximations. For concreteness, we 
assume that the point scatterers are two-level atoms each 
having a non-degenerate ground state \g) with energy E g 
and the total angular momentum J g = and an excited 
state \e) with E e = E g + Tiujq, J e = 1, and lifetime 1/T 
(h is the Planck's constant). The excited state is thus 
triply degenerate and splits in 3 sub-states with different 
projections m = — 1, 0, 1 of the angular momentum J e 
on the quantization axis z. To study light scattering in a 
random ensemble of identical atoms, we start with micro- 
scopic quantum-mechanical equations that describe the 
joint dynamics of operators corresponding to the atoms 
and to the electromagnetic field [18-20 . These equa- 
tions are solved for the field variables and equations for 
atomic operators are obtained in which coupling between 
atoms is described by the so-called "Green's matrix" G. 
It is essentially built up of Green's functions of Maxwell 
equations, describing propagation of light from one atom 
to another. G is a 3N x 3N random matrix of which a 
particular realization is determined by the ensemble of 
random positions {r^} of N atoms in 3D Euclidean space 
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Here d e i9i is the dipole moment operator of the transition 
\g) —> \e) of the atom z, = — r^, and fco = ^o/c, 
with c the vacuum speed of light. The superscripts \i and 
v denote projections of vectors on axes of the reference 
frame. 

Any excitation of the ensemble of N atoms coupled 
through the electromagnetic field can be expanded over 
eigenvectors i\) n of the matrix G. The real and imaginary 
parts of its eigenvalues A n yield the frequencies u n = 
ujq — (To/2)ReA n and decay rates T n /2 = (To/2)ImA n of 
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FIG. 1. Density of eigenvalues of the random Green's matrix. Grayscale density plots of the probability density p(A) for 
A's corresponding to long-lived states (ImA < 1). Dashed lines show the border of the eigenvalue domain following from the 
diffusion theory and the spiral branches along which eigenvalues corresponding to subradiant states are concentrated [201 121] . 
Panels (a) and (b) correspond to a low density at which the majority of eigenvalues are contained within the boundary imposed 
by the diffusion theory. Panels (c) and (d) correspond to a high density, for which states with very small decay rates ImA 
appear in the scalar model, but not in the vector one. The smallest ImA of the vector model is even larger than the prediction 
of the diffusion theory. The inset of panel (d) shows N atoms (black dots) randomly distributed in a sphere of radius R. 



the corresponding eigenstates. G is therefore the funda- 
mental object to study in order to understand the behav- 
ior of collective excitations in the atomic ensemble. We 
will be interested in the spatial localization of rj) n and 
will compare the properties of the matrix ([!]) that takes 
into account the vector character of light with those of 
its scalar approximation 

G eie . =iS eie . )^-^ (2) 

that is often used to further simplify the problem [2T|[22]. 
Indeed, it is generally believed that the vector nature 
of light does not significantly modify the main features 
of the Anderson localization problem [23 -25]. Previous 
studies of matrices similar to and ([2| supported this 
point of view [26-28 . The results that we present below 
show, however, that it is wrong for light scattering in an 



ensemble of point scatterers. 

We first analyze the density of eigenvalues A of the 
Green's matrices ([!]) and Q with a particular atten- 
tion paid to the part of the spectrum corresponding to 
long-lived states with ImA < 1. Random realizations of 
Green's matrices ([I]) and Q were generated by randomly 
choosing N points in a sphere of radius R and volume V 
[see the inset of Fig. [TJd)] . Their eigenvalues A n and 
eigenvectors i/? n = {VW;} obeying Gi/) n = A n i/? n were 
computed using Matlab or Fortran codes. Averag- 
ing was performed over 120 and 75 realizations for the 
scalar model, and over 40 and 25 realizations for the vec- 
tor model, at N = 2 x 10 3 and 4 x 10 3 , respectively. In 
the scalar case, similar results were obtained for N up 
to 10 4 . The density of eigenvalues A for N = 4 x 10 3 is 
shown in Fig. [I] At low densities p = N/V, the results 
obtained for scalar and vector models are similar, most of 
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FIG. 2. Inverse participation ratio of eigenvectors. Grayscale density plot of the average inverse participation ratio (IPR) 
as a function of the eigenvalue A of the corresponding eigenvector. Dashed lines are the same as in Fig. [I] At low density, 
subradiant states localized on pairs of closely located scatterers exist in both scalar (a) and vector (b) models. These states 
have IPR ~ 1/2. At high density, Anderson- localized states with large IPR appear in the scalar model (c), but not in the 
vector one (d). 



the eigenvalues being restricted to a region delimited by a 
line following from the diffusion theory of light scattering 
[20j[2T]. At densities exceeding p/k$ ~ 0.1, however, we 
observe that in the scalar model, a significant fraction of 
eigenvalues cross this line and acquire very small decay 
rates ImA. No such long-lived states appear in the vector 
model. 

To test the intuitive conjecture that the long-lived 
states corresponding to eigenvalues with small imaginary 
parts may be localized in space, we show in Fig. [2] maps 
of the average inverse participation ratio (IPR) for the 
same parameters as in Fig. Ill IPR n = J2iLi IVVi(ri)| 4 / 
(Y^iLi l^ni^i)] 2 ) 2 quantifies the degree of spatial local- 
ization of the eigenvector i/? n . It is of order 1/M for an 
eigenvector localized on M atoms. In the vector model, 
each VVi( r i) is a vector with 3 components and |^ n (r^)| 
should be understood as its length. As we see from Fig. 
[2j states localized on a small number of atoms exist even 



at small densities. They are typically localized on pairs 
of very closely located atoms and are due to the phe- 
nomenon of subradiance that does not require multiple 
scattering and therefore has nothing to do with Anderson 
localization [20-22 . Their eigenvalues are concentrated 
along the dashed lines that depict the evolution of the 
smallest eigenvalue of a 2 x 2 Green's matrix as the dis- 
tance between the two atoms is varied. In the scalar 
model, however, localized states of a different type ap- 
pear at densities larger than p/k^ « 0.1. These states 
have very small decay rates, in agreement with Fig. [I] 
Once again, no such localized states are seen in the vec- 
tor case. 

To convince ourselves that the localized states ap- 
pearing at large densities in the scalar model are due 
to Anderson localization, we perform the scaling anal- 
ysis [29 j. We compute the average dimensionless life- 
time of eigenstates (1/ImA) = 5uo~ l and the average 
spacing of nearest dimensionless eigenfrequencies Auj = 
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imation for £) and signals the disappearance of localiza- 
tion; the system starts to approach the effective medium 
regime. A finite- difference estimate of the scaling func- 
tion /3(g) = d\ng/d\nkoR obtained from the data of the 
main plot is shown in the inset of Fig.[3ja). As expected, 
/3(g) changes sign at g w 1, confirming Anderson tran- 
sition in the scalar model. However, none of the above 
signatures of Anderson localization is seen in Fig. |3|b) 
where we present the results for the vector model. g(k^£) 
corresponding to different N do not cross, g always re- 
mains larger than 1, and /3(g) > does not change sign, 
suggesting no localization transition. 

Our discovery of the absence of Anderson localization 
of light in a 3D random ensemble of point scatterers 
shows that completely random systems (such as, e.g., 
clouds of cold atoms) are not suitable for observation of 
this phenomenon. It also partially explains the rareness 
of its experimental realizations and indirectly supports 
the ideas about the importance of correlations in posi- 
tions of scatterers (e.g., partial order) for achieving An- 
derson localization of light [25]. However, it remains to 
be seen if such correlations, either long- or short-range, 
may induce localization. The model presented in this 
work is perfectly suitable for such a study. 

SES thanks A. Goetschy for fruitful discussions. This 
work was supported by the Federal Program for Scientific 
and Scientific-Pedagogical Personnel of Innovative Russia 
for 2009-2013 (contract No. 14.B37.21.1938). 



FIG. 3. Scaling in scalar and vector models, (a) Thouless 
number g as a function of the bare Ioffe-Regel parameter ko£ 
for the scalar model at frequency uj — ujq + To/2 and N — 
2 x 10 3 (circles) and 4 x 10 3 (squares). The curves cross at g « 
ko£ « 1 and then again at ko£ <C 1 and g « 1. Localization 
transitions take place at these points, as confirmed by the 
analysis of the scaling function /3(g) shown in the inset that 
changes sign at g ~ 1 (circles), (b) The same for the vector 
model. Solid lines in the insets are guides for the eye. 



(ReA n — ReA n _i) for eigenvalues A n in a strip of unit 
width around Re A = —1 where, according to Figs, [I] 
and [2] the localization effects are important in the scalar 
model. The localization transition for light at frequency 
uj = ujo + To/2 is expected to take place when the Thou- 
less number g = duo/Auo becomes of order unity [29| f3Q] . 
In Fig. [3] we show g as a function of the bare Ioffe-Regel 
parameter ko£, with the on-resonance mean free path £ 
calculated in the independent-scattering approximation 
(ISA) [15]: £ = kl/Aitp for the scalar and £ = k%/6np 
for the vector model, respectively. In the scalar case, 
the curves g(k^£) corresponding to different N cross at 
g « 1, ko£ « 1, as expected from the Thouless and 
Ioffe-Regel criteria of localization [2 . A second cross- 
ing takes place at much smaller ko£ (corresponding to a 
very large density p at which ISA is not a good approx- 
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